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Abstract 

Using a weighted i/ s -contraction mapping argument based on the macro-micro de- 
composition of Liu and Yu, we give an elementary proof of existence, with sharp rates of 
decay and distance from the Chapman-Enskog approximation, of small-amplitude shock 
profiles of the Boltzmann equation with hard-sphere potential, recovering and slightly 
sharpening results obtained by Caflisch and Nicolaenko using different techniques. A 
key technical point in both analyses is that the linearized collision operator L is negative 
definite on its range, not only in the standard square-root Maxwellian weighted norm 
for which it is self-adjoint, but also in norms with nearby weights. Exploring this issue 
further, we show that L is negative definite on its range in a much wider class of norms 
including norms with weights asymptotic nearly to a full Maxwellian rather than its 
square root. This yields sharp localization in velocity at near-Maxwellian rate, rather 
than the square-root rate obtained in previous analyses. 
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1 Introduction 

In this paper, we study existence and structure of small-amplitude shock profiles 
(1-1) f(x,£,t) = f(x-st,0, hm f(z) = f± 

z^±oo 

of the one-dimensional Boltzman equation 

(1-2) ft + tid x f = T- 1 Q(f,f), 

x, t € R, where f(x,t,£) G R denotes the distribution of velocities £ G R 3 at point x, t, 
t > is the Knudsen number, and 

(1.3) Q(g,h) := J (g(S , )H&)-g(Z)h&))C(Sl,t-Z m )dndt* 
is the collision operator, with 

C e R 3 , e R 3 , n e 

(1.4) e / = e+(n-(e*-o)n 

el = 6- (n.(e*-0)n. 

and various collision kernels C. Our main example is the hard sphere case, for which 

(1.5) C(fi,0 = 

See, e.g., |G1] for further details. 

Note that Q is in this case not symmetric. Other standard examples we have in mind 
are associated with the class of hard cutoff potentials defined by Grad [G] . as considered in 
[CN\ By small-amplitude, we mean that the density 



p{x,t) := (l) f {x,t) := / f{x,t,0<% 

is confined within an eo-neighborhood of some fixed reference density po > for all x, t, for 
Eq > sufficiently small, where, throughout our analysis, we have fixed 

t = 1. 

Substituting (jl.ip into (|1.2p . we seek, equivalently, stationary solutions of the traveling-wave 
equation 

(1-6) (Zi-s)d x f = Q(f,f). 

By frame-indifference, we may without loss of generality take s = 0. 
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Recall [G j IG1 [ [KMN| ICNl ILYj that the set of collision invariants (ip), that is linear forms 
such that 



t/>(OQ(jg,g)(Od£ = o, 

JR 3 

is spanned by 

(1.7) Rf:={*) f = y M/(0/(« e M 5 , *(0 = (Ui,6,£3,^| 2 ) T . 

Associated with these invariants are the macroscopic fluid-dynamical variables 

(1.8) u := Rf =: (p, pvi, pv 2 , pv 3 , pE) T , 

where p is density, v = (vi,V2,vs) is velocity, E = e + \\v\ 2 is total energy density, and e is 
internal energy density. Here, we are assuming that f(x, t, •) is confined to a space EI to be 
specified later such that the integral converges for / el. 

Taking moments of (jl.2p and applying definition (|1 .8j) . we find that the fluid variables 
obey the one-dimensional Euler equations 

Pt + d x (pvi) = 

(1.9) (pv) t + d x (vipv +pei) = 

(pE) t + d x { Vl (pE + p)) = 0, 

ei = (1,0, 0) T the first standard basis element, where the new variable p = p(f), denoting 
pressure, depends in general on higher, non-fluid-dynamical moments of /. 

The set of equilibrium states Q(f,f) = are exactly (see, e.g., (Glj) the Maxwellians 

4e/3 



(1-10) M u (0 = / , 

^(4vre/3) 3 

Making the equilibrium assumption / = M u , we obtain a closed system of equations for 
the fluid-dynamical variables consisting of the one-dimensional Euler equations (jl.9p with 
pressure p = p(p, E) given by the monatomic ideal gas equation of state 

(1.11) P =(2/3)pE. 

This corresponds to the zeroth-order approximation obtained by formal Chapman-Enskog 
expansion about a Maxwellian state [Gj IKMN] . where the expansion can be taken equiv- 
alently in powers of r, or, as pointed out in [L, MaZl], in powers of k, where k is the 
frequency in x, t of perturbations. In the present context, it is the latter derivation that is 
relevant, since (as we shall see better in a moment) we seek slowly varying solutions near a 
constant, Maxwellian, state. 

The next-, and presumably more accurate, first-order Chapman-Enskog approximation 
yields the one-dimensional Navier-Stokes equations 

Pt + d x {pu) = 

(1.12) (pv) t + d x (vxpv +p) = (p,v x ) x 

(pE) t + d x {vi(pE + p)) = (pviv x ) x + (kT x ) x , 
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where temperature T is related to internal energy by e = ^RT, R the universal gas constant, 
and 

(1.13) n = n(T) > and k = k(T) > 

are coefficients of viscosity and heat conduction. In the hard sphere case, these may be com- 
puted explicitly as /x(T) = {RT) 1 / 2 /i(l/ 7?), k(T) = (RT) 1 / 2 k(1/R) (Chapman's formulae). 
For derivations, see, e.g., [KMNj . Section 3. 

By (ll.9p , the fluid-dynamical variables associated with a traveling wave (jl.lj) must satisfy 

-sd x p + d x (pvi) = 

(1.14) -sd x (pv) + d x {vipv + pe{) = 

-sd x (pE) + d x (v 1 {pE + p)) = 0, 

hence, integrating from x = — oo to x = +oo, the Rankine-Hugoniot conditions 

(1.15) = [pvt], s[pv] = [vipv+pet], s[pE] = [vi(pE + p)], 

where [h] := h(f+) — h(f-) denotes change in h across the shock. 

Noting that endstates f± of (jl.ip by (jl.6p necessarily satisfy Q(/ ; /)± = 0, we find that 
they are Maxwellians f± = M u± , and so the associated pressures p± = p(f±) are given by 
the ideal gas formula (jl.lip . recovering the standard fact that endstates of a Boltzmann 
shock (jl.ip are Maxwellians with fluid-dynamical variables corresponding to fluid-dynamical 
shock waves of the Euler equations with monatomic ideal gas equation of state [Gj ICN) . 

This gives a rigorous if straightforward connection between Boltzmann shocks and their 
zeroth order Chapman-Enskog approximation. The following, main result of this paper 
gives a rigorous connection to the first-order Chapman-Enskog approximation given by the 
Navier-Stokes equations (j!.12[) in the limit as shock amplitude goes to zero. 

Recall [Gi] , for an ideal-gas equation of state (jl.lip under assumptions (|1.13p , that for 
each pair of end-states u± satisfying the Rankine-Hugoniot conditions (|1.15p , the Navier- 
Stokes equations (| 1 . 12[) admit a unique up to translation smooth traveling-wave solution 

u(x,t) = u NS (x - st), lim u NS (z)=u±, 

z— >±oo 

or Navier-Stokes shock. Moreover, denoting shock amplitude by e := \u+ — U-\, we have 
for e > sufficiently small the asymptotic description |Pe| 

(1.16) \d k x {u NS -u ± )\<C k e k+l e- e ^ x \ x % 0, C k , d k > 0, all k > 0. 

Up to this point in the discussion, we have made essentially no assumption on the 
nature of the collision kernel C(f2,£). For the analysis of exact profiles, we require specific 
properties of C. For simplicity of exposition, we specialize hereafter to the hard-sphere case 
(|1.5p . As discussed in Section \TT\ the arguments extend to a more general class of kernels 
including the hard cutoff potentials of Grad [G] . Then, our main result is as follows. 
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Theorem 1.1. In the hard-sphere case (jl.5p . for any given fluid- dynamical reference state 
uq and 7] > 0, there exist e$ > 0, 5 k > 0, and C k > such that for \u + — uq\ < e$ and 
e = \u + — U-\ < £q, the standing-wave equation (|1.6p /las a solution f satisfying for all 
k>0 

\d k x (u-u NS )(x)\<C k e k+2 e- 5 ^, 
(1-17) |^(/-/u NS )(x,0| < C^+V^^M^) 1 ^, 

|^(/"-/±)(x,6| < C^e-^l-lM^^) 1 ^, 

where u := i?/ is £/ie associated fluid- dynamical profile. Moreover, up to translation, this 
solution is unique among functions satisfying for < k < 2, c k sufficiently small, the weaker 
estimate 

(1-18) |^(/-/«™)(z,0| < ^e-W'l^lOs. 

Existence of small-amplitude Boltzmann profiles was established some time ago in [CN| 
for the full class of hard cutoff potentials, viewing them as bifurcations from the constant 
Maxwellian solution / = M u _ , with the somewhat weaker existence result 

\d k x (f- fu N3 )(x,0\ < C k e k+2 e- 5 ^-^M U0 (OK 
< < 1, but also the somewhat stronger result of uniqueness among solutions satisfying 

(1.19) \f-M Uo (x,0\ < Cee- s » e M-*M P M^Z)* 

for C > bounded and e > sufficiently small. For the hard sphere potential, positivity 
of profiles, and the improved estimate (11. 18ft were shown by Liu and Yu |LY| by a "macro- 
micro decomposition" method in which fluid (macroscopic, or equilibrium) and transient 
(microscopic) effects are separated and estimated by different techniques. This was used 
in |LY| to establish time-evolutionary stability of profiles with respect to perturbations of 
zero fluid-dynamical mass, J u(x)dx = 0, and thus, assuming the existence result of \CN$ . 
to establish positivity of Boltzmann profiles by the positive maximum principle for the 
Boltzmann equation (jl.2p together with convergence to the Boltzmann profile of its own 
Maxwellian approximation: by definition, a perturbation of zero relative mass in fluid- 
dynamical variables. 

The purpose of the present paper is to obtain existence from first principles by an ele- 
mentary argument in the spirit of |LY| , based on approximate Chapman-Enskog expansion 
combined with Kawashima type energy estimates [K] (the macro-micro decomposition of 
the reference), but carried out for the stationary (traveling- wave) rather than the time- 
evolutionary equations, and estimating the finite-dimensional fluid part using sharp ODE 
estimates in place of the sophisticated energy estimates of |LYj F1 In this latter part, we 
are much aided by the more favorable properties of the stationary fluid equations, a rather 

See also [Gol IHuZ] in the fluid-dynamical case. 
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standard boundary value ODE system, as compared to the time-evolutionary equations, a 
hyperbolic-parabolic system of PDE. This in a sense completes the analysis of jl.YI , pro- 
viding by a common set of techniques both existence (through the present argument) and 
(through the argument of |LY| ) positivity. At the same time it gives a truly elementary 
proof of existence of Boltzmann profiles. 

For similar results in the general finite-dimensional relaxation case, sec [MeZl, MTZ . 
The key new technical observations needed for the infinite-dimensional case are a way 
of choosing Kawashima compensators of finite rank (see Remark I4.3p . and the fact that 
the linearized collision operator remains negative definite on its range not only in norms 
of square-root Maxwellian weight where it is self-adjoint, but also in norms with nearby 
weights; this allows coordinatization with respect to a single global Maxwellian, avoiding 
unbounded commutators associated with a changing local Maxwellian frame. 

In passing, we obtain also the new result of sharp localization in velocity at near- 
Maxwellian rate (jl.17)) . which comes from improved estimates on the linearized collision 
operator independent of the basic argument. A key technical point in all three analyses- 
|CN] . |LY] . and the present one- is that the linearized collision operator L is negative 
definite on its range, not only in the standard square-root Maxwellian weighted norm, but 
also in norms with nearby weights. Exploring this issue further, we show that L is negative 
definite on its range in a much wider class of norms including norms with weights asymptotic 
nearly to a full Maxwellian rather than its square root. This observation, of interest in its 
own right, yields through the same existence argument sharp localization in velocity at 
near-Maxwellian rate, rather than the square-root rate obtained in previous analyses. 

Finally, we note that stability of small-amplitude Boltzmann shocks has been shown in 
[LY] with respect to small H s perturbations with zero mass in fluid variables. It would be 
very interesting to continue along the same lines to obtain a complete nonlinear stability 
result as in [SX] or [MaZl[ IMaZ2| . with respect to general, not necessarily zero mass, 
perturbations. 

2 The nonlinear collision operator 

We begin by a careful study of the collision operator. Related results may be found, for 
example, in [Cj [GPS] . 

2.1 Splitting of the collision operator 

In view of definition (11. 3h , we split 



(2.2) Q + (g,h)(0 = / Qn(g,h)dQ, Q a (g,h) = / gtf)h&)C(n,t - ?)d& 




Q(g,h) = Q+(g,h)-Q-(g,h) 



into gain and loss parts [G] [Gl], where, for Q defined as in (II. 4p 
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and 



(2.3) Q-(g,h) = g(OMO, MO = J c(n,e. - Oh(Z*)dt*dn. 



2.2 Estimates for Q_ 

In the hard sphere case (jl.5 



(2.4) v h {£) = J C(n,^-t)h(^)d^dn = c J \£-n)\ h(rf)dr,. 

Here and elsewhere, denote (£) := (1 + l^] 2 ) 1 / 2 following standard convention. 

Lemma 2.1. In the hard-sphere case (|1.5p . for h > with (Qh G L , is positive, 
continuous and 

(2-5) (0<^(0<<0|MCi- 
Proof. Evidently, 

^(0<A|e|+B, A = \\h\\ L i, B=\\\ V \h\\ L t. 
This implies the upper bound. Next, 

mo > a\h\ - b 

which implies the lower bound for |£| > (B+l) /2A. For £ bounded, the integral is continuous 
and bounded from below. □ 

2.3 Estimates for Q + 

Consider the Maxwellians 

(2.6) W .(0 = e-^l 2 

and the weighted I? spaces HP = w s L 2 (M 3 ) with norm 

(2-7) ||/|U = (/e 2 ^ 2 |/(6| 2 ^)" 

We use similar definitions for spaces Hf [resp. 1H| ] of functions of one variable £i [resp. 
two variables (£2, £3)]- 



Proposition 2.2. In the hard-sphere case (|1.5p . for any s > 0, 

(2-8) \\(0'^Q+(9,h)\\ w < IMUKO^IU 

(2-9) ||Q+(^fc)|| H , < IMUKOMU- 



We first estimate Qq (see definition (|2.2p ) for a fixed fisS 2 . 
Lemma 2.3. In the hard-sphere case (|1.5f) . for s > 0, 

( 2 - 10 ) \\(Z)-*Qn(<j,h)\\ m < C s \\g\\ ms \\(0^h\\ ms 

and 

( 2 - n ) \\Qn(9,h)\\ m < C s \\g\\ w \\(0h\\ w 

Proof. Choose the orthonormal basis such that 

n = (1,0,0). 

In this case 

£ = (£*i,£2)£3)j = (£i,£*2,£*3)) C(^, £* — £) = l£i ~~ 

and has the more explicit form 

(2-12) Ql(9,h)(0 = [ g(vi^2,^)h(^m,V3)\^-Vl\dr] 



In particular 

(2-13) Qi(9,h)(0=Ig(Si>&,&)Jh(Si) 
with 

(2.14) 6,6) = J\Zi-m.\§(rn,t2,&)dm 

(2.15) A(6) = y Kii,m^m)dmdm- 

Note that J(£i) = (£i)~^e-"l^l 3 j($i) with 

m) = (Zi)* J e- s(v > W3 \l + g + vl + vir l *H(Zi,V2,V3)dV2d V3 
with ||iI|| L 2 = 11(0 ^||h». Thus, 

\\3\\l* < \\(0^h\\ w . 

Therefore, 
with 

r If, _ le-slml 5 
G(0 = y H^rnWrj^&drn, k(^m) ■= — ^ 



and \\g\\ L 2 < \\g\\w°- Note that 



sup ||A;(£i, -)|| L 2 = C < oo. 
61 



Thus 

G(ei,6,6) 2 <C 2 0(6,6), \\4 L1 < \\g\\l 2 

It follows that 

\\q\\ L 2 < C\\j\\ L2 \\g\\ L2 . 

The proof of the second estimate is similar □ 

Integrating over Q on the unit sphere, implies Proposition ^. 21 Combining the estimates 
for Q_ and Q + implies: 

Corollary 2.4. In the hard-sphere case (ll.5p . for all s, there is C s such that 

(2.I6) \\{^Q{g,h)\\ ns < C s ||<OMUI<ONU 

2.4 Further estimates 

Proposition 2.5. In the hard sphere case (|1.5|) . suppose that < s < s' , and g £ M s and 

h G HP' . Then 

(2-17) ||Q+(5,/i)|| HS < C S)S /||5|| HS ||/t|| Hs ,. 

Proof. Introduce G = e s ^ 2 \g(£)\ and H = e s '^ 2 \h(()\. Then, when Q = (1,0, 0) T , 

91 :=e 2s ^\Q n (g,hm 

satisfies 

gi (0<e (s - s ' )lal Y G (^'^,e3)ff(ei,%,%)iCi-%|e"^ 12 ^ 

By the Cauchy-Schwarz inequality, 

kl(6| 2 <(ei) 2 e 2(W)la|2 | |G(^,6,6)| 2 \H(^, m , m )\ 2 d V . 
where we have used that 

J Iei-r7i| 2 e- 2 ^l 2 dr/<(ei) 2 
Since s' > s, (^i) 2 e 2 ^ s ~ s '^^' 2 is bounded and 

\\ qi \\ L2 < \\G\\ L2 \\H\\ L2 . 

Integrating over f2 implies the result. □ 
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Proposition 2.6. In the hard-sphere case () 1 . 5 1) . suppose that < s < s' , and g 6 I s and 
heW. Then 

(2.18) ||Q + (5,h)|| H . <C«,«'||^|| H .'||^|| H .. 

Proof. Introduce G = e s '\^ 2 \g{£,)\ and H = e s ^ 2 \h(£)\. Then 

qn:=e 2 ^ 2 \Q n (g,h)(0\ 

satisfies 

qn(0 < e ('-* , NKI a -«-n) 3 ) J *si(Z,tiM - r,) • n\e- s ^ 2 drj 
where for all fi£S 2 : 

ll < ML 2 (IR3 xR 3) < || G\\ L 2 \\H\\ L 2 

Integrate over Q and use the Cauchy-Schwarz inequality to get 

\e 2s ^ 2 Q + (g,h)(0\ 2 < | J qn(H)dn\ 2 <aJ \^,r,)\ 2 d V dn 



with 

A 

Thus 



: = J e 2 ( s - s ')(isiM^)i 2 )| ( e _ n) ■ n\ 2 e-^\ 2 d v dn. 

A<1 + | e 2 ( s - s ')(iei 2 -(C-n)| 2 )| e . fi |2 d0 . 



To compute the integral (for large £), we can choose coordinates such that £ = (0,0, t), 

■ 5 tt, 2?rJ, 



t > 0, and parametrize the sphere with angular coordinates 9 G [0, 27r] and <p 6 [— h7t, i7r] 



so that Q = (cos p cos 0, cosy? sin 9, sin 99) In this case the integral becomes 

2vr J** e 2{s ~ s "> t2cos2 n 2 sin 2 ip cos if dtp 



which is smaller than 

1 . 



2vr f'^ e 2{ - s - s '^ 2cos2 ^2t 2 s\npcosp dip 
Jo 

2tt C e 2 ^ t2u t 2 du = 2vr f e 2 ^*'^ < 



Therefore j4 is uniformly bounded and 

| |e 2s KI 2 Q + (<7,>0(6| 2 ^< / \^n)\ 2 dr}dm^<\\Gf L2 \\Hf L ,. 



□ 



Corollary 2.7. in the hard-sphere case fl 1 . 5 [> . suppose that < s < s' . and f £ M s . Then 
for a£ff f/ie mappings f 1— > Q-(a, /), / 1— > Q+(a, /) and / Q+(f, a) are bounded from 
M s to M s ; with norm controlled by a constant times ||a|| Hs /. 
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Remark 2.8. The estimates above were proved for convenience for the Gaussian weights 
u) = &€\ and uj s . They immediately extend to any Maxwellian weight M u and M 3 . 

3 The linearized collision operator 

We next study the linearized collision operator about a Maxwellian or nearby velocity 
distribution. Fix a reference state u. The associated Maxwellian M R is denoted by M_. 
For s s]0, 1], let M s denote the space of functions / on M 3 such that 

(3-1) = / mO~ 2 V(C)\ 2 dC < +oo. 

Note that M e HP for all s < 1. The space Hs plays a particular role as it will be clear 
below. 

Proposition 3.1. The quadratic mapping f i— > Q(f,f) is continuous from (^)~^M S to 
keri?n ((^alP), where R is the operator (|1.7|) defining the thermodynamical variables. . 

Proof. The action from (£) _ 2]H S to (£)5IHP is a consequence of Corollary 12.41 and Re- 
mark 12. 81 That the image is contained in ker R follows from the known properties of the 
collision operator: 

(3.2) V/eff: RQ(fJ) = 0. 

□ 

Given a function a, the linearized collision operator at a is 

(3.3) L a g = Q(a,g) +Q(g,a). 

In particular, we consider first the linearized operator at a = M_: 

(3.4) Lg = Q'M9 = Q(M,g)+Q(g,M). 
Corollary 12.71 implies the following result: 

Lemma 3.2. For all s £ L is a bounded linear operator from (£) _ 2]HP to keri?n 

{£)m s and from {^^M 3 to keriinlEP. 

3.1 Symmetry and coercivity on ElK 

Let V = ker Rf]M.2 and let U denote the orthogonal complement of V in H. It has dimension 
5. Noting that 

(3-5) (x,f) L2 = (xK,f) R 
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we see that U is spanned by the functions vpjM_. An orthogonal basis is 

(3-6) MO = Xj(OM.(0, 3 = 0,..., 4, 

with 

Xo(0 = l. Xf(0 = ^7=J* for j = 1,2,3, 
(3-7) : ^ 

We denote by Pu and Py the orthogonal projection from Hs to U and V respectively. 



In the language of |LYj . U is the macroscopic part of / and V is the microscopic part. 
Note that Uc(0^ffi and Ocff for all s < 1. Therefore 

Lemma 3.3. For s E [±,1[, P v maps IP onto V s :=YnH s = ker J RnM s and (£) 5HP onto 

vn ((O^ihp). 

Remark 3.4. The projections Pu and Py do not commute with the operator of multipli- 
cation by £i. They are not orthogonal in HP for s > 5, but still produce a continuous 
decomposition H s = U © V s . 

Proposition 3.5. L is (formally) self adjoint and nonpositive in H2 and definite negative 
on V. More precisely, 

i) for all f and g in f £ (£) _1 Hl2, 

(3-8) {Lf, 9 )^ = {f,Lg)^. 

ii) there is S > such that for all f G (^) _1 H: 

(3-9) S\\(0^vf\\ 2 mh <-Re(Lf,f) mi . 

Notes on the proof. This is a classical result in the theory of Boltzmann equation in the hard 
sphere case and more generally in the case of hard cut off potentials (see e.g. [C j lG ] IG1 [ ICN] ) 
1. The analysis of section [2] splits L into L = + K, with 

Kg(0 = - J fci(£,£*)0(&R* + / fe(^fl.e*)j(0^<fii + / h(^,n,^) g (Qd^dn 

with 

fe2(e.n.e.)=M(Oc(n,e,-0 
fc 3 (C,^^)=M(Oc(^C*-6 
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Using the conservations 



(3.10) e+& = iei 2 + i^i 2 = ief + ie:i 2 , = 

which imply that 

(3.H) M(«) = M(OM(e), 

one shows that 

{M- 1 K j g,h) L2 (g,Mr 1 K J h) L2 

for j = 1, 2, 3, implying the symmetry of L in H2. 

2. One can also argue as follows. By Boltzmann's H -theorem, 



/ 



Q(/,/)logM<0 



for all / withe enough decay at infinity. Hence, Taylor expanding about the Maxwellian M, 
a minimizer of J Q(f, /) log fd£,, we obtain symmetry and nonnegativity of the Hessian, 

(Lh)h , 

giving nonnegativity of L on H2 and also formal self-adjointness. 

3. It is known that is K compact in EI 2 and that kerL = U (this can be proved using 
the formulas above). By self-adjointness of L on HI 2, to establish strict negativity on V, 
it is sufficient to establish a spectral gap between the eigenvalue zero and the essential 
spectrum of L. But, this follows from Weyl's Lemma by comparison of L = —v + K with 
the multiplication operator by —v{£) < — Co < 0. □ 



In a more explicit form, the inequality (|3.9h reads 

(3-12) 5 J (OMtfrWCOl 2 ^ < -Re J M(t;)- l Lf(0f(m- 

We also point out the following properties which are freely used below and which follow 
from the symmetry of L in HI 2 : 

(3.13) L = ¥ y L = LF Y , F V L = LF V = 0. 

3.2 Coercivity on M s 

With A > to be determined later on, introduce the equivalent norms: 
(3-14) \\f\\l--=\\f\W+ X \\f\\U- 
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(3.16) 



Proposition 3.6. For \ < s < 1, the operator L is continuous from (£)~2W S to ((£) 2HP)nV 
and from (£)~ 1 H. S to V s , and formally coercive on V s for the norm (|3.14p . More precisely, 
there are A > and 5 > such that for all f 6 (^) _1 V S : 

(3-15) «5||(0ip v /||| s <-Re(L/,/) gs 
Proof. We want to prove that 

^^(M^^ + am^)- 1 )!/^)! 2 ^ 

< _ Re y (M(0- 2s + AM^)- 1 )^/^)/^)^, 
Following the analysis of Section [21 

L = -v ld + K 

where fo(0 ~ (0 an d ^ is bounded from HP to HP, since the Maxwellian M 6 BP' for 
s < s' < 1. Hence there is <5i > such that 

*i / (0M(0- 2s i/(0i 2 ^< | M(e)- 2s L/(o/(e)^ + c | M(o- 2s i/(oi 2 rfe 

Moreover, by (|3.12p . there is 5o such that 

So J (omo~ l \f(o\ 2 <% < j mo^Lfionodt;- 

Hence: 

r (0{tiM(0~ 2s + A^oM^)" 1 )!/^)! 2 ^ 



< -Re (L/,/) gs +C y M(0- 2s |/(e)| 2 ^- 
We choose A such that for all £ 

CM(£)~ 2s < i(e)(5iM(e)- 2s + AJoM^r 1 ) 

implying the inequality (|3.16j) , □ 

Remark 3.7. Included in the bound (13. 15H is the observation that both the first-order 
Chapman-Enskog approximation t/jvs and the entire hierarchy of higher-order Chapman- 
Enskog correctors lie in HP, any < s < 1, something that is not immediately obvious. 
Indeed, looking closely at the inversion of L a , we see that they in fact decay at successively 
higher polynomial multipliers of the full Maxwellian rate. 



15 



3.3 Comparison 

We consider the linearized operator L a g = Q(a,g) + Q(g,a) at a, not necessarily nonnega- 
tive, close to M . 

Proposition 3.8. For s G [|,1[ and a G (^"llfF, L a is bounded from (£) - tlrP to (£)iIEF. 
Moreover, there are constants 5 > 0, C > 0, A > and £o > such that for all f G (()~2i s 

(3.i7) ||(e)-^a/||i3<c||(e)^p v /|| fis + c £ ||(e)^Pu/|| fr , 

(3-18) | (AJ,/) gs | < C||(e)3Pv/|||. +Ce 2 \\(0^ " u ' 2 

and 

(3.19) *||<0W||1- <-Re(L a /,/) as + C £ 2 ||(O^I 
with 

(3.20) e = e(a):=\\(Oka-M)hs- 

In (|3.19p the H s scalar product has to be understood as the integral 

(3.21) (lj, /) g . = y (M(o- 2s + wto-^Lammodt, 

which is well defined since / G (£)~2HP and L a f G (()^H S . 

Proof. That L a is bounded from (£)~2H S to (£) 2ff for all s G [i, 1[, follows directly from 
Section [2j Moreover, 

L a f-Lf = Q(a-M,f) + Q(f,a-M) 

and 

\\(0~kLaf-Lf)\\^<Ce\\(Oh\\ns- 
Since Lf = LPyf, this implies (|3.17p . Since L a f and Lf belong to V and thanks to the 
definition of the modified scalar product HP, 

[L a f-Lf,f)z, = (LJ-Lf,F v f)^ = o(s\\(Oh\\ m \\(0^f\ 



With (|3.17|) . this implies (|3.18p with a new constant C. With Proposition 13.61 this implies 
(I3TT9]) . □ 

Remark 3.9. Since U is finite dimensional, one can use any norm for Fjjf in the estimates 
(13371) flgJgH and (I3TT9D above. 

Remark 3.10. We have in mind that e{a) can be taken arbitrarily small. This holds if 
a = M u and u is close to u since , when s < 1, 

\\(0HMu-M)\\ 2 ms = J(0\M u -M\ 2 M- 2s dH 

as |u — u| — > by Lebesgue's dominated convergence theorem. 
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4 Abstract formulation 



We now rephrase the problem in a general framework, for the square-root Maxwellian norm 
H = H2 in which we carry out the main analysis. We treat general weights in Section 
110.21 by a bootstrap argument. Taking the shock speed equal to by frame-indifference, 
we consider (II. 6h as the abstract standing-wave ODE 

(4.1) AU' = Q(U,U). 
with 

(4.2) Af(0 = 6/(0 

independent of U (semilinearity of the Boltzmann equation), and Q as in (|1.3j) . (11. 5h . 

4.1 Bounds on the transport operator 

The collision operator has been studied acting in spaces HP associated to our reference 
Maxwellian M. We have the following evident facts regarding the transport operator A. 

Proposition 4.1. For s £ [i,l[, the operator A is bounded from {^sff to (£)~2lHP and 
(formally) self adjoint in HP as well as in HP . 

4.2 Kawashima multiplier 

We next construct a Kawashima compensator as in [K, MeZl, MTZ], but taking special 
care that the operator remain bounded in this infinite-dimensional setting. 

Proposition 4.2. There are C , 5 > 0, A > and there is a finite rank operator K 6 
Jz^HI -1 ,!! 1 ) such that such K is skew symmetric in HP and satisfies 

(4.3) Re (KA - L) > 5(£}Id. 
meaning that 

(4-4) c||(0'/|||.<Re((A'A-L a )/,/) g .<C||(03/|||.. 

Proof, a ) We first check that the genuine coupling condition is satisfied, i.e. that there is no 
eigenvector of A in ker L = U. Indeed, using the basis <pj of U given in (13. 6h . an eigenvector 
of A with eigenvalue r in U is a linear combination ^ such that the polynomial 

4 

(&-r)X>iXj-(0 
j=o 

is identically zero. Equating to zero the term of degree 3 implies that 04 = 0. Equating to 
zero the coefficient of the terms of degree 2 implies that ctj = for j = 1, 2, 3, and finally 
ao = 0. Thus the property is satisfied. 
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b) We look for K as 
(4.5) K = e(K u +K 12 + K 21 ) 

with > a parameter to be chosen and 

K 12 = A l2 := ¥ V AF Y = A* 2l , K 21 = -K* l2 = -¥ Y A¥ V := -A 2l , 

and 

K n = Pu#iiPu = -^ii- 

Here * means the adjoint with respect to the scalar product in H. We have used Proposi- 
tion ED 

Thus, with An '■- 



(4.6) Re FvKAPv = ± [K n ,A U } + AX\A 2X , 

The condition a) means that An (restricted to U) has no eigenvector in ker A 2 i = ker A\\A 2 i , 
with A 2 iA 2 i symmetric positive semidefinite and An symmetric. Since dimU is finite (equal 
to 5), this implies by the standard, finite-dimensional construction of Kawashima et al [K] 
that one can choose Kn such that Re FuKAFu is definite positive on U: there is ci > 
such that 

(4.7) (Re P v KAF v f,P v f) m > ci||Pu/|£. 
Moreover, since dimU is finite, there is another c\ > such that 

2 



(4.8) (Re P v KAP v f,¥ v f) n > c t \ r |lh , . 

Thus, using Proposition 13.91 for a = M: 



2 . II™ ,||2 

H 1 



Re [(KA-L)f,f) m >9ci\\P v f\\" ml + c 

- 0C l|/|| H i 
with 

C = ||ifiiPu^Pv|| + ||JCi2Pv^Pv|| + ||i^2iFu^Fu|| + ||K2iFu^F v || 

where the norms are taken in J£ (H ; HI -1 ). All these operators have finite rank < n and 
are bounded. Thus if 9 is small enough, this shows that Re (KA — L) is definite positive 
in the sense of ()4.4p . Using the perturbation Lemma 13.81 implies that the estimate remains 
true for a satisfying (|3.20p . □ 

Remark 4.3. The construction above, by reduction to the equlibrium manifold, is essen- 
tially different from the original proof of [K] in the finite-dimensional case, which would 
yield a symmetrizer of infinite rank. The advantage of finite rank is that we need not worry 
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about boundedness of the operator. We note that this is related to methods in the Boltz- 
mann literature in which the Kawashima compensator is replaced by estimates on a reduced 
Chapman-Enskog approximation such as the Grad 13 moments model or the Navier-Stokes 
approximation, again to avoid possible boundedness issues; see, e.g., [G j ILY| . 

See also the related construction of |GMWZ| in the case that u is scalar, for which K\\ 
may be taken equal to zero. We note that we could apply the same reduction argument 
to the reduced problem and proceed by iteration to this scalar case, thus obtaining an 
alternative proof in the finite-dimensional case as well. 



4.3 Reduction to bounded operators 



In the hard-sphere case (II. 5h . we may rescale the equations to obtain a problem involving 
only bounded operators. We have EI 1 C EI C EI -1 , bounded operators from H 1 to EI -1 and 
we work with the scalar product of H. We can multiply the equations on the left by 
if A G if (H 1 ; H _1 ) then 

A := (0~ X A G ^(H 1 ; H 1 ) 

and 

(Afj) m = (Af,f) m 

so that if A is symmetric in H, A is symmetric in H 1 . 

Equivalently, we can make the change of variable / i— > / = from H 1 to H and 

define 

Af:= (C)- l 2Af = (Cr^A(C}- l 2f- 

Then A G jSf (H; H) and A is symmetric in H if A is. 

By Corollary 12.41 the corresponding collision operators 

Q(f,f) ■= (O^QifJ) 

and 

Q(fJ) ■■= 

by Corollary El are bounded as well: Q G ^(H 1 ; H 1 ) and Q G 3§(M;M), where 3S{H;H') 
denotes the space of continuous bilinear forms from H — > H', i.e., B G 3)(H\H') if and 
only if 



(4-9) \\B{g,h)\\ H , <C s \\g\\ H \\h\ 



ii ■ 



4.4 The framework 

At this point, we have reduced to the following abstract problem, with semilinear structure 
quite similar to that treated in the finite-dimensional analysis of [MeZlj . Working in EI 
with operators A and Q and dropping tildes, we study the standing-wave ODE 

(4.10) AU' = Q(U,U), 

with U taking its values in an infinite dimensional space EI. 
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4.4.1 Assumptions on the full system 

We make the following assumptions, verified above for the Boltzmann equation in the hard- 
sphere case with A, Q replaced by A, Q. 

Assumption 4.4. (i) A is a bounded self adjoint operator in a (real) Hilbert space HI; 
(ii) There is an orthogonal splitting H = U © V with U finite dimensional 
(Hi) Q is bilinear, continuous (in sense (|4.9p ) and symmetric (in U) from Ixl toY. 

For U E H, we denote by Ljj the bounded operator V i— > 2Q(U, V), that is the differential 
of Q(U, U). We denote by Pu and Py the orthogonal projectors on U and V respectively. 
We use the notations U = u + v, with u = F\jU and v = FyU . 

Assumption 4.5. We are given a reference state U_ (in a smaller space M C HJ such that 
L = Ljj is self adjoint with kernel U and L is definite negative on V. 

Lemma 4.6. There are 5 > 0, Eq and C > such that for a E M and U E H: 

(4.11) - Re (L a U, U) R > 5\\F Y U\\ 2 m - Ce\\F v U\\ m \\F y U\\ m 
provided that 

(4.12) ||a-M|| H <e<e - 
Proof. By continuity of Q, (|4.9p . there is C such that 

(4.13) \\L a U-Lu\\ m < C||a-M|| H ||Z7|| H . 
Moreover, there is 5 > such that 

-{LU,U) m = -{LU,F Y U) m = -{LF Y U,F Y U) m > S\\¥yU\\^. 

Since 

(L a U,U) m = Re (L a U,F Y U) m 
the lemma follows. □ 

Lemma 4.7. In an ^.-neighborhood of U_, the zero set of Q is given by a smooth (indeed 
C°°) manifold M = {U : v = v*(u)} with : U — > V smooth. 

Proof. Assumption 14.51 and the Implicit Function Theorem, together with the observation 
that Q as a continuous biinear form (in sense (|4.9p ) is C°° in the Frechet sense. □ 



We further assume the Kawashima condition established in Proposition 14.31 

Assumption 4.8. There is a skew symmetric bounded operator K E Jf(M) and a constant 
7 > such that 

(4.14) Re KA - L > 7 Id. 
Using (|4.13p . this implies 

Lemma 4.9. There are 7 > and eq > such that for a E H satisfying (|4.12|> and U E FJ; 

(4.15) Re((KA-L a )U,U) m > 1 \\U\\ 2 m 
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4.4.2 Assumptions on the reduced system 

Coordinatizing U G H as 



(4.16) U 
we have 



u\ (P V U 
v \PvU 



(4-17) A=(> t 2 V Q ' ° 



a 21 A 22 y -* v?^) 

with An E JSf(U;U), Ai 2 G JSf(V;U) etc. We use the notation 
(4.18) /i(u, u) := iuM + A i2 v G U. 

Finally, the equilibria are parametrized by u: 

u 



(4.19) M w- V »,(«y 

where is the smooth mapping from a neighborhood of u to a neighborhood of v = v*(u) 
in V, as described in Lemma 14.71 

Recall from jY] that the reduced, Navier-Stokes type equations obtained by Chapman- 
Enskog expansions are 

(4.20) K(u)' = {h(u)u'Y, 
where 

(4.21) h*(u) := h(u, v*(u)) = A\\U + Ai 2 v*(u), 

(4.22) h(u) := -A 12 c*(u) 



c*(«) := d v q («, w*(u)) 

A21 + A 22 cfo*(u) - dv*(u)(Au + j 4i 2 g?v*('u) ) . 



with 
(4.23) 

Note also, by the Implicit Function Theorem, that dv*(u) = —d v q~ 1 d u q(u,v*(u)). 

An important property of the Chapman-Enskog approximation, following either by di- 
rect computation or by coordinate-independence of the physical derivation, is that the form 
(|4.2ip - (|4,23p of the equations is coordinate invariant, changing tensorially with respect to 
constant linear coordinate changes; moreover, the change in functions h*, 6* due to a con- 
stant linear coordinate change may be computed directly from (|4.20p using the coordinate 
change in u alone. From this we find in the Boltzmann case that (|4.20j) is equivalent through 
a constant linear coordinate change to the Navier-Stokes equations (|1.12p with monatomic 
ideal gas equation of state and viscosity and heat conduction coefficients satisfying (|1.13D . 
We make the following assumptions on the reduced system, verified for the Navier-Stokes 
equations (hence satisfied for the Boltzmann equation) in [MaZ3j . 
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Assumption 4.10. (i) There exists s{u) symmetric positive definite such that sdh* is 
symmetric and sb* is symmetric positive semidefinite. 

(ii) There is no eigenvector of dh* in ker . 

(Hi) The matrix b*{u) has constant left kernel. 

(iv) For all values of u, ker ir*dh*(u) n kerb*(u) = {0}, where tt*(u) is the zero eigen- 
projection associated with b*(u). 

Finally, we assume that the classical theory of weak shocks can be applied to (|4.20p . 
requiring that the flux /* have a genuinely nonlinear eigenvalue near 0. 

Assumption 4.11. In a neighborhood U* of a given base state uq, dh* has a simple eigen- 
value a near zero, with a(«o) = 0, and such that the associated hyperbolic characteristic 
field is genuinely nonlinear, i.e., after a choice of orientation, Va • r(uo) < 0, where r 
denotes the eigendirection associated with a. 

Remark 4.12. As discussed in [Yj, Assumptions I4.10( i)-(ii) hold in great generality. 
Assumptions I4.107 iii)-(iv) must be checked in individual cases. 

4.5 The basic estimate 

With these preparations, we can establish existence by an argument almost identical to 
that used in [MeZl] to treat the finite-dimensional case: indeed, somewhat simpler. The 
single difference is that in carrying out the basic symmetric energy estimates controlling 
microscopic variables we do not attempt to exactly symmetrize L a at each x value as 
was done in [ MeZlj . but only use the fact that each L a is approximately symmetric by 
construction. This is important in the infinite-dimensional case, since exact symmetrization 
can (and does in the Boltzmann case) introduce unbounded commutator terms that wreck 
the argument. To isolate this important technical point, we carry out the key estimate here, 
before describing the rest of the argument. 
We consider the equation 

(4.24) A8 x U-L a U = F 
with a = a{x) satisfying 

(4.25) \\d k Mx) -M)\\ m < C k e k+l e-^\ 
We assume that 

(4.26) F V F = eh + d x f. 

Lemma 4.13. There is a constant C such that for e sufficiently small, one has 

(4.27) \\U'\\ L 2 + \\¥ Y U\\ L2 < C(\\f\\ H 2 + \\h\\ H i + \\g\\ H i + e\\¥ v U\\ L2 ). 
Here, the norms L 2 , H 1 etc denote the norms in L 2 (M;EI), ff 1 (lR;IHI) etc. 
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Proof. Introduce the symmetrizer 

(4.28) S = d 2 x + d x o K - Aid. 
One has 

Re dl(Ad x - L a ) = -Re d x oL a od x - Re d x o L dxa 
Re d x o K(Ad x — L a ) = d x o Re o 9 Z — Re d x o If L a 
Re (A^ - L a ) = -Re L a 

where ReT = ^(T + T*) and the adjoint is taken in L 2 (M; H). We have used that [d x ,L a ] = 
Lg xa by linearity of L with respect to a. Thus 

Re S o (Ad x - L a ) =d x o (Re AK" - L a ) od x - ARe L a 
- Re d x o L 9a:a - Re <9 X o iiTL a . 

Therefore, for U G F 2 (M), (|4.11j) . <|4.15|> and the continuity of K and Q imply that 

Re (SF,U) L 2 > 7 \\d x U\\l 2 + X[6\\F v U\\l 2 - Ce\\F v U\\ L 2\\F v VU\\ L 2) 
- Cll^allL-ll^llL^I^^IlL 2 - C\\d x U\\ L 2\\L a U\\ L 2. 

We note that 

(4.29) L a U = LF V U + (L a - L)U 
Therefore, 

(4.30) \\L*U\\hh % PvU\\hh + 4^ U \\hh 
Taking A large enough and using (|4.25p yields 

\\U'\\ 2 L2 + \\¥ V U\\ 2 L 2 < Re (SF,U) L 2+e\\F v U\\ L 2(\\FvU\\ L 2 + \\U'\\v). 
In the opposite direction, 

Re (SF,U) L 2 ^[IS^H^dlS^lUa + 11^1111^11^) 

+ A(e||/ l ||£2||Pu«|U» + \\f\\ L 4Wd x U\\ L 2 

+ \\F Y F\\ L 2\\F Y U\\ L 2). 

The estimate ()4.27p follows provided that e is small enough. 

This proves the lemma under the additional assumption that U G H 2 . When U G H 1 , 
the estimate follows using Friedrichs mollifiers. □ 



5 Basic L 2 result 

We now describe a simpler version of our main result, carried out in the L? norm H. For 
clarity of exposition, we carry out the entire argument in this more transparent context, 
indicating afterward in Section [10] how to extend to the general (pointwise, higher weight) 
norms described in Theorem ll.il 
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5.1 Chapman— Enskog approximation 

Integrating the first equation of (|4.1U|) and noticing that the end states (u±,v±) must be 
equilibria and thus satisfy v± = v*(u±), we obtain 

/r A 11 u + Ai2V = f*(u±), 

(5.1) , 

A 2 \u + A22V =q(u,v). 

Because / is linear, the first equation reads 

(5.2) f*(u)+A 12 (v-v*( y u)) = f4u±). 

The idea of Chapman-Enskog approximation is that v — v*(u) is small (compared to the 
fluctuations u — u±). Taylor expanding the second equation, we obtain 

(A21 + A 2 2dv*(u))u + A 2 2(v - v*(u))' = d v q(u,v*(u))(v - v*(u)) 

+ 0(\v-v*(u)\ 2 ), 

or inverting d v q 

v — v*(u) = d v q~ 1 (u,v*(u)){A2i + A2 2 dv*(u))u' 
+ 0(\v - v*(u)\ 2 ) + 0(\(v - 17,(U))'|). 

The derivative of (|5.2p implies that 

(A n u + Ai 2 dv*(u))u' = 0{\{v - v*(u))'\). 

Therefore, f|5.3|) can be replaced by 

(5.4) v - v,(u) = c,(u)u' + 0(\v - v,{u)\ 2 ) + 0{\{v - 

where c* is defined at f|4.23[) . Substituting in (|5.2|) . we thus obtain the approximate viscous 
profile ODE 

(5.5) 6„(u)u' = /.(u) - /*(u±) + 0(|u - v*(u)\ 2 ) + 0(|(« - u,(it))'|), 

where 6* is as defined in (|4.22|) . 

Motivated by (|5.3p - (|5.5p . we define an approximate solution (uns,vns) of <|5.1|) by 
choosing TL]yg as a solution of 

(5.6) b*(v,Ns)u' NS = f*(u NS ) - f*{u±), 
and vns as the first approximation given by (|5.3|) 

(5.7) u^S - v*(nArs) = c*(u NS )u' NS . 

Small amplitude shock profiles solutions of (|5.6|) are constructed using the center mani- 
fold analysis of |Pe| under conditions (i)-(iv) of Assumption 14. 10( see discussion in [MaZ5] . 
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Proposition 5.1 ( |MaZ5j ). Under Assumptions ^. 10 and \4.1i , in a neighborhood of(uo,uo) 
in R ra x M. n , there is a smooth manifold S of dimension n passing through (uq,uq), such 
that for (u_,u + ) G S with amplitude e := \u + — u_| > sufficiently small, and direction 
(u+ — U-)/e sufficiently close to r(uo), the zero speed shock profile equation (|5.6p has a 
unique (up to translation) solution Uns inU*. The shock profile is necessarily of Lax type: 
i.e., with dimensions of the unstable subspace of dh*(u-) and the stable subspace of dh*(u + ) 
summing to one plus the dimension of u, that is n + 1. 

Moreover, there is 9 > and for all k there is C k independent of (u__,u + ) and e, such 
that 

(5.8) \d k x {u NS - u ± )\ < C k e k+1 e~ e ^, x^O. 

We denote by S + the set of (u-,u + ) G S with amplitude e := \u + — u_| > sufficiently 
small and direction (u+ — u_)/e sufficiently close to r(u$) such that the profile u~ns exists. 
Given (u-,u + ) G S + with associated profile u~ns, we define vns by (|5.T[) and 

(5.9) U NS := (u NS ,v NS ). 

It is an approximate solution of (|5.1[) in the following sense: 
Corollary 5.2. For (u_,u + ) G 5 +; 

(5.10) A n u NS + A 12 v NS - f*(u±) = 
and 

7^ := + ^22^5 - q(u~NS,VNs) 

satisfies 

(5.11) |d^(x)| < C k £ k+3 e- eeM , x^O 
where C k is independent of (n_, u+) and £ = \u+ — u~\. 

Proof. Given the choice of vns, the first equation is a rewriting of the profile equation (|5.6p . 
Next, note that 

vns ~ v*(u NS ) = 0(\u' NS \), (v NS - v*(u NS ))' = 0(\u'xs\) + 0(l^7vsl 2 )' 

where here O(-) denote smooth functions of uns an d its derivatives, which vanish as indi- 
cated. With similar notations, the Taylor expansion of q and the definition of vns show 
that 

K v =0(\v NS - v*(u NS )\ 2 ) + 0(\(v NS - v*(u NS ))'\) 
+ dv*(u NS (A n + A 12 dv*(u N s))u' NS - 

Moreover, 

(A u + Ai 2 dv*(u N s))u' NS = {f*(uNs)) = {b*(uNs)u NS )) 

= 0(\u' NS \ 2 )+0(\u^ s \). 

This implies that 

K v = 0(\u' NS \ 2 ) + 0(K s \). 
satisfies the estimates stated in (|5.1ip . □ 
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Remark 5.3. One may check that if we did not include the correction from equilibrium 
on the righthand side of (|5.7p . taking instead the simpler prescription v^s = v*(uns) as 
in [LYj . then the residual error that would result in (|5.10|) would be too large for our later 
iteration scheme to close. This is a crucial difference between our analysis and the analysis of 
|LYj . The prescription Ujys corresponds to the first-order Chapman-Enskog approximation 
in both variables, u and v together. 



5.2 Basic L 2 result 

We are now ready to state the basic L 2 version of our main result. Define a base state 
Uq = (uo,v*(uq)) and a neighborhood U = x V. 



Proposition 5.4. Let Assumptions (SS), (GC), and \4-10 hold on the neighborhood U of 



Uq, with Q G C°° . Then, there are Eq > and 5 > such that for (u_,u + ) & S+ with 
amplitude e := \u + — u—\ < Eq, the standing-wave equation (j4.10jl has a solution U in IA, 
with associated Lax-type equilibrium shock (u—,u+), satisfying for all k: 

\d k x (U -U NS )\<C k e k+2 e- 5e \ x \, 
(5.12) \8 k (u - u±)\ < C k E k+l e~ S£ \ x \ x ^ 0, 

\d k x {v -v*{u) I < C k E k+2 e- 5£ W, 

where Uns = (unSjVns) is the approximating Chapman-Enskog profile defined in (15. 9j) . 
and C k is independent of e. Moreover, up to translation, this solution is unique within a 
ball of radius ce about Uns in norm \\ ■ \\ L 2 +E^ 1 \\d x -\\ L 2 +E~ 2 \\d 2 ■ \\ L 2 , for c > sufficiently 
small. (For comparison, Uns ~ U± is order e 1 / 2 in this norm, by (|1.17p (ii)-(iii).) 



6 Outline of the proof 

We describe now the main steps in the proof of Proposition 15.41 exactly following the finite- 
dimensional analysis of [MeZlj . 



6.1 Nonlinear perturbation equations 

Defining the perturbation variable U := U — Uns, and expanding about Uns, we obtain 
from (|5.ip the nonlinear perturbation equations 

(6.1) A n u + A 12 v = 

(6.2) A 2 iv! + A 22 v' - dq{U NS )U = -1Z V + N{U) 

where the remainder N(U) is a smooth function of Uns an d U, vanishing at second order 
at U = 0: 

(6.3) N(U)=M(U NS ,U) = 0(\U\ 2 ). 



26 



We push the reduction a little further, using that 



(6.4) M := dq(u NS ,v NS ) - dq(u NS ,v*(u NS )) = 0(\v NS - v*(u NS )\)- 

Therefore the equation reads 



(6.5) 





-K v + MU + N(U) 



where 

(6.6) Q 2 i = d u q(u N s,v*(u N s)), Q22 = d v q(u NS ,v*(u NS ))- 
Differentiating the first line, it implies that 

(6.7) LlU:=AU'-dQ(u N3 M*Ns))U= (_ nv + MU + N{U) 

The linearized operator Ad x — dQ(U) about an exact solution U of the profile equations 
has kernel U 1 , by translation invariance, so is not invertible. Thus, the linear operators 
L% and C% are not expected to be invertible, and we shall see later that they are not. 
Nonetheless, one can check that L% is surjective in Sobolev spaces and define a right inverse 
£*(£*)^ = or solution operator {£%)* of the equation 



(6.8) ClU = F 



as recorded by Proposition 16.21 below. Note that L% is not surjective because the first 
equation requires a zero mass condition on the source term. This is why we solve the 
integrated equation (|6.5p and not ()6.7j) . 

To define the partial inverse (£*)' , we specify one solution of (|6.8p by adding the co- 
dimension one internal condition: 

(6.9) 4 • u(0) = 
where l e is a certain unit vector to be specified below. 

Remark 6.1. There is a large flexibility in the choice of l E . Conditions like (|6.9p are known 
to fix the indeterminacy in the resolution of the linearized profile equation from (|5.6|) and it 
remains well adapted in the present context, see section [8] below. A possible choice, would 
be to choose t e independent of e and parallel to the left eigenvector of dh*(uo) for the 
eigenvalue (see Assumption 14. 1 1]) . which, together with the asymptotics of Proposition 
15.11 gives 

(6.10) 4 • U' NS (0) ~ £ 2 + 0. 
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6.2 Fixed-point iteration scheme 

The coefficients and the error term 7Z V are smooth functions of uns and its derivative, thus 
behave like smooth functions of ex. Thus, it is natural to solve the equations in spaces 
which reflect this scaling. We do not introduce explicitly the change of variables x = ex, 
but introduce norms which correspond to the usual H s norms in the x variable : 

(6.11) 11/11^=^11/11^ +s-Hd x f\\ L2 + ... + £ l-||^/|| i2 . 

We also introduce weighted spaces and norms, which encounter for the exponential decay 
of the source and solution: introduce the notations. 

(6.12) < x >:= (x 2 + 1) 1/2 

For 8 > (sufficiently small), we denote by s the space of functions / such that e Se<x> f G 
H s equipped with the norm 

(6.13) ||/||^ =e^c~ k \\e 5£<X> d k J\\L>- 

k<s 

Note that for S < 1, this norm is equivalent, with constants independent of e and 5, to the 
norm 

\\ n Se<x> f\\ 

\\ e jWhi- 

Proposition 6.2. Under the assumptions of Theorem \5.4\ there are constants C, eo > 
and 5q > and for all e G]0, eo], there is a unit vector i £ such that for e G]0, eo]; <5 S [0, So], 
f G s , g G ff^ the operator equations (|6.8j) (|6.9p has a unique solution U G H 2 s , 
denoted by U = (£%yF, which satisfies 

( 6 - 14 ) IIO^L* < Ce- l {\\f\\ H , + \\g\\ H *). 

Moreover, for s > 3, there is a constant C s such that for e G]0, eo] and f G , 
g G H s e 5 the solution U = (C^F G H s e 5 and 

(6.15) HOC^Hh. < C e - 1 (||/|| iys + 1 + )+C s \\(jO^F\\ He - 1 . 

The proof of this proposition comprises most of the work of the paper. Once it is 
established, existence follows by a straightforward application of the Contraction-Mapping 
Theorem. Defining 

( 6 - 16 ) T --=^{-n v + Mu + Niu)^ 

we reduce (]6.7p to the fixed-point equation 
(6.17) TU := U. 
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6.3 Proof of the basic result 

Proof of Theorem \5.4\ The profile uns exists if e is small enough. The estimates 
imply that 

(6.18) \\u NS - u±\\ H « iS < C s e 

with C s independent of e and 8, provided that 8 < 6/2. Similarly, (|5.1ip implies that 

(6.19) \\n v \\Hi s <C s e 3 , 
and (|6,4p implies that 

(6.20) \\M\\ H s s <C s e 2 . 

Moreover, with the choice of norms (|6,lip . the Sobolev inequality reads 

(6.21) < C\\u\\ H i < C\\u\\ H i 

£ £ ,<5 

with C independent of e. Moreover, for smooth functions there are nonlinear estimates 

(6.22) PWIIffj < C(\\u\\ L °o) \\u\\ Hi . 
which also extend to weighted spaces, for 8 < 1: 

(6.23) Mu)\\ Hls < C(\\u\\ Loo ) \\u\\ Hls . 

In particular, this implies that for s > 1, 8 < min{l,0/2} and e small enough: 

\\MU\\b U < C(\\M\\ Hl jU\\ Hls + ||M||^J|C/||^) 
<e 2 (C\\U\\ H s s + C s \\U\\ Hls ) 

where the first constant C is independent of s. Similarly, 

(6-25) \\N(U)\\h 1s < C{\\U\\ l ~)\\U\\ h1s \\U\\h U . 

Combining these estimates, we find that 

WTUWh' < e-\C s e* + Ce 2 \\U\\ H s + C s e 2 \\U\\ H , + C\\U\\ H i \\U\\ H 'X 

that is 

(6.26) < C s £ 2 + C( £ + ^ 1 H^II^,)II^II^,, + C s e\\U\\ HL 

provided that £ < 8 < min{l,#/2} and ||C/||loo < 1. 
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Consider first the case s = 2. Then, T maps the ball B e 5 = {\\U\\ H 2 < s 1+ i} to itself, 

e,S 

if e < Ei where e± > is small enough. Similarly, 

(6.27) \\TU-TV\\ H 2 s <Ce-\e 2 + \\U\\ m + \\V\\ H 2)\\U -V\\ Hh , 

provided that ||{7||z,°° < 1 and || V||i,oo < 1, from which we readily find that, for e > 
sufficiently small, T is contractive on B £ g, whence, by the Contraction-Mapping Theorem, 
there exists a unique solution U £ of (|6.17p in B £: $ for e sufficiently small. 
Moreover, from the contraction property 

\\U £ - T(0)\\ H? = \\T(U £ ) - T(0)\\ m < c\\U £ \\ H? , 

with c < 1, we obtain as usual that ||C/ £ ' 5 || 0-2 < C||7~(0)|| u-2 , whence 

£,5 £,5 

(6.28) \\U £ \\ H 2 <Ce 2 . 

by (|6.26l) . In particular, e eS ^U £ = 0(e 2 ) in H 2 and by the Sobolev embedding 

(6.29) \\e e5{x) U £ \\ L <~ =0{e 2 ), \\e £S{x) d x U £ \\ L ~ = 0(e 3 ). 

For s > 3, the estimates (I6.26D show that for e < £\ independent of s, the iterates T ra (0) 
are bounded in H^ s , and similarly that T ra (0) — T(0) = 0(e 2 ) in H^ s , implying that the 
limit U belongs to iT? s with norm 0{e 2 ). Together with the Sobolev inequality (I6.21h . this 
implies the pointwise estimates ()1.17j) . 

Finally, the assertion about uniqueness follows by uniqueness in £> C£j( 5 under the addi- 
tional phase condition (|6.9p for the choice 5 = and c > sufficiently small (noting by 
our argument that also B CEt s is mapped to itself for e sufficiently small, for any c > 0), 
together with the observation that phase condition (|6.9[) may be achieved for any solution 
U = U NS + JJ with 

ll^'IU- < ce 2 « U' NS {0) ~ £ 2 
by translation in x, yielding U a (x) := U(x + a) = Uns(x) + U a (x) with 

U a (x) := U NS (x + a) - U NS {x) + U(x + a) 

so that d a (£ £ ■ u a (0)) = £ £ ■ (u' NS (a) + u'(a)) ~ l e ■ u' NS (0) and so (by the Implicit Function 
Theorem applied to h(a) := £~ 2 {i e • u a ), together with l e ■ uq = o(e) and the assumed 
property that t a • u' NS (Q) ~ e 2 coming from our choice of £ £ ; see (16.10j) . Remark l6.ip the 
inner product t a ■ u' NS (0) may be set to zero by appropriate choice of a = o(e~ 1 ) leaving U a 
in the same o(e) neighborhood, by the computation U a — Uq ~ d a U ■ a ~ o(s~ 1 )s 2 . □ 

It remains to prove existence of the linearized solution operator and the linearized bounds 
(|6.15p . which tasks will be the work of most of the rest of the paper. We concentrate first 
on estimates, and prove the existence next, using a viscosity method combined with (the 
single new step in treating the infinite-dimensional case) discretization in velocity. 
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7 Internal and high frequency estimates 



We begin by establishing a priori estimates on solutions of the equation (|6,8p This will 
be done in two stages. In the first stage, carried out in this section, we establish energy 
estimates showing that "microscopic" , or "internal" , variables consisting of v and derivatives 
of (u, v) are controlled by and small with respect to the "macroscopic" , or "fluid" variable, 
u. As discussed in Section [4.51 this is the main new aspect in the infinite-dimensional case. 

In the second stage, carried out in Section [8l we estimate the macroscopic variable u by 
Chapman-Enskog approximation combined with finite-dimensional ODE techniques such 
as have been used in the study of fluid-dynamical shocks |MaZ4l IMaZ5j IZTj |Z2"1 IGMWZj . 
exactly as in the finite-dimensional analysis of [MeZl| . 



7.1 The basic H 1 estimate 

We consider the equation 

(7.1) ClU : 

and its differentiated form: 

7' 



A u u + A 12 v *) = (f 

A 2 iu' + A 22 v' - dq(u NS ,v*(u NS ))U J \g 



(7.2) AU' - dQ(u NS , (u NS ))U = C 



The internal variables are U' = (u',v') and v where 

(7.3) v:=v+pu, p = d v q~ 1 d u q(uNS,v*(uNs)) = -dv*(v,Ns) 

is the linearization about (ujvs, vns) of the key variable v — arising in the Chapman- 
Enskog expansion of Section [5TT1 Noting that pu = at the reference point U by Assumption 
31 we have the important fact that 



(7.4) IM|h = OOON|h 

on the set of U we consider (e 2 close to Uns, so e close to U), so that v and v are nearly 
equivalent. 

Proposition 7.1. Under the assumptions of Theorem \5.4\ for there are constants C, Eq > 
anddo > such that for < e < Sq andO < 5 < So, f G H^ 5 , g G s andU = (u,v) G H^ s 
of (|7.ip satisfies 

( 7 - 5 ) H^ILa +IHLa <C{\\(f,f'J",g,g')\\ L 2 +e\\u\\ L2 ). 

Proof. For 5 = 0, the result follows by Lemma 14.131 together with (|7.4|) . 
For 5 > small, consider U w = e e5 ^ U. Then, U w satisfies 



(7.6) L%U % 



r 



9 W 1 



31 



with f w = e eS ^f and g w = e e5 ^g + ed(x)' \A 21 u w + A 22 v w ). We note that, 
\\U'\\ L 2 < \\(U™)'\\ Li +e\\ir\\ I z„ \\v\\ L 2 < \\v w \\ L >, 

\\r,(n',(n , '\\Li<\\(fj , ,n\\L>, 

£ £,d 

\\g w ,(9 w y\\Li<\\(9,g')\\ L ^ + ^\\(u,u')\\ L2 

£ £,d £,d 

We use the estimate (|7.5p with 5 = for f/ 1 ", and the Proposition follows provided that 8 
is small enough. □ 

7.2 Higher order estimates 

Proposition 7.2. There are constants C, Eq > 0, 5q > and for all k > 2, there is 
such that < e < Eq, 5 < 5q, U E H^ s , f E H^ 1 and g E H s £ $ satisfying (|7.2j) satisfies: 

\\d k x U'\y +\\d k x v\\ L 2 < C\\d*(f,f'j",g,g')\\ L 2 

+ e k C k (\\U'\\ H k-i + e[|5||„*- 1 + e\\u\\ L 2 J 

s,S e,6 

Proof. Differentiating (|7.ip k times, yields 



(78) Ad x U k - dQ(u NS , v*(u NS ))d x U - , „ fc 



where 

fc-i 



d£ 5 + r k 



u x V22 ^ " 
Z=0 

Here we have used that dq(uNS,v*(uNs)U = Q 22 v. The -ff 1 estimate yields 
\\d k x U'\\ L , + \\d k x v + V d k x u\\ L , < C(\\d k (fj'j",g,g')\\ L2 

+e||^«|| L 2 +||<V fc || L 2 +||r fe || i2 ), 

£,0 £,d £,0 

for < < s, with ro = when k = 0. Since Q is a function of u^s, its /c — Z-th derivative 
is 0(e fc ~ /+1 ) when k — I > 0. Therefore: 

Rr fe || i2 +||r fc || i2 < C k £ k (\\v'\\ Hk ^ +e\\v\\ L 2 ). 

Similarly, for = 1 

ll^fcll^ < \\d x v + pd x u\\ L 2 + Ce 2 \\u\\ L 2 

e,o £,o e,o 

and for k > 2: 

\\d k v k \\ L 2 < ||^ + p^u|| i2 +C jfc (£ fe ||n / ||^- 2 +e fe+1 ||n|| i 2 ). 

e,o e,o e,S 

□ 
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8 Linearized Chapman— Enskog estimate 

8.1 The approximate equations 

It remains only to estimate ||u||r2 in order to close the estimates and establish (|7.5p . 
To this end, we work with the first equation in (|7.ip and estimate it by comparison with 
the Chapman-Enskog approximation (see the computations Section 15. lj) , exactly as in the 
finite-dimensional case |MeZl| . 
From the second equation 

A 2 iu' + A 2 2v' - g = d u qu + d v qv = d v qv, 

where we use the notations v of Proposition 17.11 we find 

(8.1) v = d v q~ l f(i 2 i + A 22 d v dv* (u NS ))u f + A 22 v' - g 
Introducing v in the first equation, yields 

(An + A 12 dv*(u N s))u + A 12 v = /, 

thus 

(A n + A 12 dv*(u NS ))u' = f - A 12 v' - d 2 v*(u NS )(u' NS ,u). 
Therefore, (|8.1|) can be modified to 

(8.2) v = c*(uns)u' + r 
with 

r = d^qiuNSiV^UNs)) (a 22 (v)' - g 

+ dv*(u NS )(f - A 12 v' - d 2 v*(u NS )(u NS ,u))y 
This implies that u satisfies the linearized profile equation 

(8.3) b*u' — dh*u = A\ 2 r — f 
where = K(uns) and dh* := dh*(uNs) = An + Ai 2 dv*(v,Ns)- 

8.2 L 2 estimates and proof of the main estimates 

The following estimate was established in [MeZl| using standard finite-dimensional ODE 
techniques; for completeness, we recall the proof here as well, in Section [8.31 below. 

Proposition 8.1 ( |MeZl| ). The operator b*d x — dh* has a right inverse (b*d x — dh*)^ 
satisfying 



(8.4) IKhdx-dh^hW^KCe- 1 
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uniquely specified by the property that the solution u = (b*d x — dh*)^h satisfies 
(8.5) 4 • u(0) = 0. 

for a certain unit vector l e . 

Taking this proposition for granted, we finish the proof of the main estimates in Propo- 
sition 16.21 



Proposition 8.2. There are constants C, Eq > and 5q > such that for e G]0, Eq], 
5 G [0,do], / G H^s, g G H^ s and U G H^ s satisfying JO} and 

(8-6) II^^C^dl/llH^ + NI^ J- 

£,o t >° £,o 

Proof. Going back now to (|8.3|) . it satisfies 

- d\u = 0{\v'\ + \g\ + |/'| + e 2 |u|) - /, 
If in addition u satisfies the condition (|8.5|) then 

(8-7) \\u\y < Cs-\\\v'\\l* + ll(/)/'>5)IL2 +e 2 ||n|| L 2 ). 

£,0 £,0 £,0 £,d ' 

By Proposition 17.11 and Proposition 17.21 for k = 1, we have 
( 8 - 8 ) ll^'ll^ + < C(ll (/,/', /"i 5, 5') || L 2 + e ||«|| z2 ). 

£,5 £,<5 £ '" £.(5 

H^'IL* +||«1L 2 < 

(8 9) 

' C7(||(/',r 5 r )ff ',/)|| i2 +,||C/'|| L2 +e 2 \\u\\ L2 ). 

Combining these estimates, this implies 

ll ll-t^.5 Ml eS ll sS ll ll-^ e( 5 / 

<C(e|| (/,/',/", ff^Ol^a +e 2 ||n|| L2 ). 
Substituting in (|8.7p , yields 

£||n|| i2 < C(||(/,/', 5 )|| L2 +e||(/,/',/",fl,</)lltfi +s 2 || u || i2 ). 

£,0 v £,0 £,0 £,0 7 

Hence for e small, 

(8-10) e\\u\\ Lh <C{\\(fJ\g)\\ Lh +e\\(fJ\f\g,g')\\ Hl J. 

Plugging this estimate in (|8.8[) 



(8.11) \\U'\\ L2 +\\v\\ l2 +e\\u\\ L2 <C\\(f,f',f",g,g')\\ H i +). 

£,5 £,5 £,(5 
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Hence, with (18.91). one has 



\\U"\\ L 2 +\\v'\\ L2 < 

1.12) 



C{\\(f'J",r,g',g")\\ L z+e\\(f,f,f",g,g')\\ H i). 



e,5' 

Therefore, 

(8.13) \\U'\\ H1 +\\v\\ L 2 +44l* <C\\f,f',f",g,g'\\ H i 

E,5 £,S £,S £,S 

The left hand side dominates 

II ^11 Hi + £ \\ U '\\l 2 = £ \\ U '\\h 2 

£,5 £,5 £,5 

and the right hand side is smaller than or equal to ||/||„ 2 + <7 H i ■ The estimate (I8.6P 

£,§ £,S 

follows. □ 
Knowing a bound for ||u||r2 , Proposition 17.21 immediately implies 

£,S 

Proposition 8.3. There are constants C , £q > and 5q > and for s > 3 there is a 
constant C s such that for e G]0, £q\, 5 G [0, 5q], f G H s e + S , g G s and U G s satisfying 
()6.8p and ()6.9p . one has 

(8.14) ||tf|| ^Ce^iWfW^+r + WgW^ )+C s \\U\\ HS - 1 . 

£,6 £,& £,6 £,6 

Remark 8.4. The estimate of Proposition 18.11 may be recognized as somewhat similar 
to the estimates of Goodman |Go] obtained by energy methods in the time-evolutionary 
case, the same ones used by Liu and Yu [LY] to control the macroscopic variable u. More 
precisely, the argument is a simplified version of the one used by Plaza and Zumbrun [PZ] 
to show time-evolutionary stability of general small-amplitude waves. 

8.3 Proof of Proposition 18.11 

By Assumption l4.10T i). we may assume that there are linear coordinates u = (ui,u 2 ) G 
R ni x W 12 and h = (hi,h 2 )e K ni x W 12 , with n 2 = rank b*(u) such that 

(8-15) K(u) = ( 

\b 21 (u) b 22 [u) 

and b 22 (u) is uniformly invertible on U*. Introducing the new variable 
(8.16) u 2 = u 2 + V Ul , V = {b 22 )- l b 21 {u NS ), 

the equation b*u' — dh*u = h has the form: 

a ll u\ + a 12 u 2 = hi, 

(^•l^) 722 -/ -21 -22- t 

b u 2 — a u\ — a u 2 = h 2 
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where 



„ , Id 0\ r {0 



Assumption 14. lOT ii) implies that the left upper corner block o 11 is uniformly invertible. 
Solving the first equation for u\ , we obtain the reduced nondegenerate ordinary differential 
equation 

bf u' 2 + a 21 (a lv r l a 12 u 2 - a 22 u 2 = h 2 + a 21 (a n r% 

or 

(8.18) bu' 2 -au 2 = h = 0{\h 1 \ + \h 2 \). 

Note that det dh* = det d 11 det a by standard block determinant identities, so that 
detd ~ det dh* by Assumption 14. lO( ii) . Moreover, as established in [MaZ4], by Assumption 
14.111 and the construction of the profile uns we find that m := (6) _1 d has the following 
properties: 

i) with m± denoting the end points values of m, there is 9 > such that for all k : 

(8.19) \d k x (m(x)-m ± )\<e k+1 e-^- 

ii) m(x) has a single simple eigenvalue of order e, dented by efi(x), and there is c > 
such that for all x and e the other eigenvalues A satisfy |Re A| > c; 

iii) the end point values [i± of fi satisfy 

(8.20) > a n+ < -a 

for some a > independent of e. 

In the strictly parabolic case det 6* 7^ 0, this follows by a lemma of Majda and Pego 
[MP] . 

At this point, we have reduced to the case 

(8.21) u' 2 -m{x)u 2 = 0{\h l \ + \h 2 \), 

with m having the properties listed above. The important feature is that m' = 0(e 2 ) « e, 
the spectral gap between stable, unstable, and e-order subspaces of m. The conditions 
above imply that there is a matrix uj such that 

p := w _1 mo) = blockdiag{p + , efi,p— }, 

where the spectrum of p± lies in ±Re A > c. Moreover, uj and p satisfies estimates similar 
to (I8.19p . The change of variables u 2 = ujz reduces ()8.2ip to 

(8.22) z - pz = uj- 1 uj'z + 0(\hi\ + \h 2 \). 
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The equations (z + )' — p + z + = h + and (z~)' — p~z~ = h~ either by standard linear 
theory [He] or by symmetrizer estimates as in [GMWZJ , admit unique solutions in weighted 
L 2 spaces, satisfying 

||e 6|a: U ± || i 2 < CHe^l^H^, 

provided that 6 remains small, typically 6 < |Re p |. 

The equation z' — efizo = ho may be converted by the change of coordinates x — ► x := ex 

to 

(8.23) dxz - fl(x)z = h (x) = £ -1 /to(x/e), 
where Zq(x) = zo(x/e) and jl(x) := n(x/e). By (|8.19l) 

\fi(x)-fi ± \ < Ce- e ^ 

with n± satisfying (|8.20p . This equation is underdetermined with index one, reflecting the 
translation-invariance of the underlying equations. However, the operator d$ — fl has a 
bounded I? right inverse (d$ — 1 , as may be seen by adjoining an additional artificial 
constraint 

(8.24) 5b(0) = 

fixing the phase. This can be seen by solving explicitly the equation or applying the gap 
lemma of [McZ2] to reduce the problem to two constant-coefficient equations on x ^ 0, with 
boundary conditions at z = 0. We obtain as a result that 

He 51 * 1 z \\ L i<C\\e 5 Wh \\ L 2 

if 5 < min{a, 6}, which yields by rescaling the estimate 

H„ed|x| „ II ^ 1 1| „£<J|a;| u II 

II e ' ' zq\\ ^2 < u e || e 'nolli 2 

Together with the (better) previous estimates, this gives existence and uniqueness for 
the equation 

z — pz = h, zo(0) = 

with the estimate ||e e<5 ' :r ''z|| i 2 < Ce -1 ]^ 515 ^^!!^. Because uj~ x uj' = 0(e 2 ), this implies that 
for e small enough, the equation (18.221) with zq(0) = has a unique solution. Tracing back 
to the original variables u, the condition zq(0) = translates into a condition of the form 
t e ■ u(0) = 0. Therefore, the equation b*u' — df*u = h has a unique solution such u that 
l £ ■ u(0) = 0, which satisfies 

\\e e5 ^u\\ L 2 < Ce^We^hWv 
for 5 and e small enough, finishing the proof of Proposition 18.11 
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9 Existence for the linearized problem 

The desired estimates ()6.14j) and (|6.15p are given by Propositions 18.21 and 18.31 It remains 
to prove existence for the linearized problem with phase condition u(0) ■ r(e) = 0. This we 
carry out using a vanishing viscosity argument. 

Fixing e, consider in place of C%U = F the family of modified equations 

(9.1) CVW-ClU-vK) =F:= ((], £ £ -u(0) = 0. 



Differentiating the first equation yields 

(9.2) ALT' - dQ(x)U -U»=( f \ 4 ■ u(0) = 0. 
where dQ{x) denotes here the matrix dQ{u^ g ,v*(u~n s)) ■ 

9.1 Uniform estimates 

We first prove uniform a-priori estimates. We denote by S" the Schwartz space and for 
S > 0, by ^ e s the space of functions u such that e eS ^u G S", with (x) = \/l + x 2 as in 

dsn. 

Proposition 9.1. There are constants Eq > 0, 5q > and rjo > 0, and for all s > 2 

a constant C s , such that for e G]0,£o]; S £ [0, So], rj G]0,r/o], and U and F in y e $(R), 
satisfying (I9.1h 

(9.3) \\U\\ H . ^Cse-WlfWjjs+i + WgW^ ). 

e,5 ■'■ " - s 



Proof. The argument of Proposition 17.11 goes through essentially unchanged, with new r\ 
terms providing additional favorable higher-derivative terms sufficient to absorb new higher- 
derivative errors coming from the Kawashima part. 

Thus we are led to equations of the form (|7.2p with the additional term —rjU" in the left 
hand side. Using the symmetrizer S (|4.28j) . one gains ^||^"||^2 +A||C/'||^ 2 in the minorization 
of Re (SF, U) and loses commutator terms which are dominated by 

^ll^ , ||Lo (||[/'||2 2 + || [ /|| i2 ||[/'|| i2 ) + ^H^ll^dl^ll^ + || [711^)11 C/"|| i2) 
which can be absorbed by the left hand side yielding uniform estimates 

(9.4) Vv\\U"\\ L 2 + \\U'\\ L 2 + \\v\\ L 2 < C{\\f\\ H 2 + \\h\\ m + \\~g\\ m +e\\u\\ L 2). 
Going back to (|9.2p , this implies uniform estimates of the form 

(9.5) Vv\\U"y +\\U'\\ L2 +\\v\\ L2 <C{\\(fJ',f",g,g')\\ L2 +e\\u\\ L2 ). 

e,S s,S £ >° £,S 

for 5 = 0, and next for 5 E [0, 5q] with Sq > small, as in the proof of Proposition 17.11 
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When commuting derivatives to the equation, the additional term rjd 2 brings no new 
term and the proof of Proposition 17.21 can be repeated without changes, yielding estimates 
of the form 

s/v\\D k x U^ Lh +\\d k x U'\\ Lh + \\d k x v\\ Lh 

(9.6) <C\\d k x (f,f'J",g,g')\\ Lh 

+ e k C k (\\U'\\ H k-i +e\\v\\ H k-i +s\\u\\ L 2 ). 

s,5 e,6 £ >° 

Next, applying the Chapman-Enskog argument of Section [8] to the viscous system, we 
obtain in place of (|8.3|) the equation 

(9.7) hu' - dKu = f + 0(\v'\ + \g\ + |/'|)+ e 2 0(\u\) + r]0{\u'\ + \U"\), 

where the final r\ term coming from artificial viscosity is treated as a source. One applies 
Proposition 18. II to estimate e||it||£2 by the l? e ^-norm of the right hand side, and continuing 

as in the proof of Proposition 18.21 the estimate f)8. 13j) is now replaced by 

Vv\\u"'\\ Lh + \\u'\\ Hls + \\v\\ L 2 g +4 u Li s 

(9 ' 8) < C(\\f,f',f",g,g'\\ Hl + V (W'\\l* + W"\\ L *)). 

Therefore, for r/ small, the new 0(rj) terms can be absorbed, and (|9.3p for s = 2 follows as 
before. The higher order estimates follow from (19. 6ft . □ 



9.2 Existence 

We now prove existence and uniqueness for (|9.ip . First, recast the the problem as a first- 
order system 

(9.9) W - MA = T 

with 




and 

/ An A12 \ 

(9.10) A := rf l rjl \ . 

\n- l A 2 iA ll - Q 21 r ) - 1 A 2 xA 12 - Q 22 A 22 ) 

Next, consider this as a transmission problem or a doubled boundary value problem on 
x ^ 0, with boundary condtitions given by the n + 2r matching conditions U(0~) = U(0 + ) 
at x = together with the phase condition i e ■ u(0) = 0, that is n + 2r + 1 conditions in all: 

(9.11) W(0 _ ) =W(0+), 4-«(o) = o. 
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Note that the operator-valued coefficient matrix A converges exponentially to its end- 
states at ±00, by exponential convergence of Uns and boundedness of A, Q. 

Lemma 9.2. There is 9\ > such that for e small enough , the limiting coefficient matrices 
A± have no eigenvalue in the strip |Re z\ < e5q. 

Proof. The proof is parallel to the proof of the estimates. Dropping the ±, suppose that ir 
is an eigenvalue of A, or equivalently that there is a constant vector U 7^ such that e lTX U 
is a solution of of equations (|9,ip Thus 

A\\u + A12V = irrju, 

(9 12) 

(it A — Q + T 2 rj)U = 0. 

In the first equation, introduce once again the variable v = v + Q22Q21U, so that the 
equations are transformed to 

A* n u + A 12 V = ITTjU, 

{it A - Q + t 2 t])U = 0. 

Denoting by K the end point values of the Kawashima multipliers associated to A and Q, 
consider the symmetrizer 

E = |t| 2 — ifK - A. 

Multiplying the second equation in (|9.13p by S and taking the real part of the scalar product 
with U yields 

|r| 2 Re (KA - QU, U) + X(QU, U) + v\t\ a {U, U) 

<C{\lmr\(\T\ 2 + X))\U\ 2 + C\T\\QU\\U\ 

+ r/(|r| 2 |Imr| 2 + |r| 3 + A|r| 2 ))|[/| 2 . 

Therefore, choosing appropriately A, for 77 and |Imr| sufficiently small, one has 

(9.14) (?7|r| 4 + \t\ 2 )\U\ 2 + \v\ 2 < C(|Imr| + e)|n| 2 

In particular, |r| must be small if Imr, e are small. 

From the equation irA2iu+ A22V— Q 21U— Q 22V +rjT 2 v = and the fact that IQ21I = 0(e) 
by (|7.4p . one deduces that 



iT(Q 2 2)~ 1 A 21 u = 0(\t\ +t}\t 



2m„~. 



V . 



Substituting in the first equation of (|9.13|) . we obtain the Chapman-Enskog approximation 

(A* u - irb*)u = 0(t]\t\ + |r| + 7/|r| 2 )|Imr|^))|u| 
where 6* denotes the end point value of the function (|4.22p . Therefore, 
(9.15) {(b^A^u - iru\ < C|Imr|i |r||«| 



40 



with arbitrarily small c > 0. We know from Assumption 14.111 that for e small, (fr*) -1 ^^ 
has a unique small eigenvalue, of order 0(e), real. Let us denote it by ejj,. Then we know 
that |^| is bounded from below, see (I8.20p . Then (19.150 implies that there is a constant C 
such that for |Imr| small enough, and thus |r| small, |ir — efi\ < C|Imr|2 |r|. Therefore, 
/Imr + efi\ < ^e\fi\ if e is small enough. 

Summing up, we have proved that if e is small enough, A has at most one eigenvalue z 
in the strip |Rez < e2|/i|, such that \z — e\i\ < k £ \fJ-\- This implies the lemma. □ 

Remark 9.3. The same reasoning can be applied to prove that A actually has a simple 
eigenvalue such that \z — £ji\ < ^e|^|. 

9.2.1 Finite-dimensional case 

We first review the case that U is finite-dimensional, recalling for completeness the analysis 
of [MeZlj . 

Proposition 9.4 ([MeZlj). There are constants Eq > 0, 5q > and t]q > such that 
for e g]0,£o]j 8 £ [0, <5o], V 6]0 3 Vo}> an d F ^ n ^e<s(K), (|9.ip admits a unique solution 

Proof. Noting that the coefficient matrix A converges exponentially to A± at ±oo, we may 
apply the conjugation lemma of [McZl] to convert the equation (|9.9[) by an asymptotically 
trivial change of coordinates U = T{x)Z to a constant-coefficient problems 

(9.16) Zi-A_Z_ = F_, Z' + -A + Z + = F + , 

on {±x > 0}, with n + 2r + 1 modified boundary conditions determined by the value of the 
transformation T at x = 0, where A± := A(±oo), and Z±(x) := Z(x) for ±x > 0. 

By standard boundary- value theory (see, e.g., [He]), to prove existence and uniqueness 
in the Schwartz space for the problem (19. 9\\ on {x < 0} and {x > 0} with transmission 
conditions (|9.1ip . it is sufficient to show that 

(i) the limiting coefficient matrices A± are hyperbolic, i.e., have no pure imaginary 
eigenvalues, 

(ii) the number of boundary conditions is equal to the number of stable (i.e., negative 
real part) eigenvalues of A + plus the number of unstable eigenvalues (i.e., positive real part) 
of A_, and 

(iii) there exists no nontrivial solution of the homogeneous equation / = 0, g = 0. 
Moreover, since the eigenvalues of A± are located in {|Rez| > 9i£, the conjugated form 

(|9.16p of the equation show that if the source term / has an exponential decay e~ e5 ^ 
at infinity, then the bounded solution also has the same exponential decay, provided that 
5 < 6\ . Therefore, the three conditions above are also sufficient to prove existence and 
uniqueness in S?^ if e and 5 are small. 



41 



Note that (i) is a consequence of Lemma 19. 2| while (iii) follows from the estimate (|9.3p 
To verify (ii), it is enough to establish the formulae 



(9-17) , , „ , , ,„ 



dimS(A ± ) = r + dim5(A*f ) 
dim^(A ± ) = r + dim^(A*± 



where A^f = dh 1r {u±) = An + A\2dv*(u±) and S(M) and U{M) denote the stable and 
unstable subspaces of a matrix M. We note that A^f = dh*(u±) are invertible, with 
dimensions of the stable subspace of A^ and the unstable subspace of A*^[ summing to 
n + 1, by Proposition 15.11 Thus, (|9.17p implies that 

dimS(A + ) + dim^(A_) = 2r + dim<S(A^) + dimU(A^) = 2r + n + 1 

as claimed. 

To establish (|9.17p . introduce the variable v = v + Q22Q21U, and the variable corre- 
sponding to v 1 scaled by a factor 772, that is w = rfzw + 77^2 Q^Qixi^Awu + A\2v). After 
this change of variables, the matrix A it conjugated to A with 




1 T / .1 I 

2 




(9.18) r/5A 



From (i), the matrix 7/ 2 A has no eigenvelue on the imaginary axis, and the number of 
eigenvalues in {Re A > 0} is independent of rj, and thus can be determined taking rj to 
infinity. The limiting matrix has r eigenvalues in {Re A > 0}, r eigenvalues in {Re A < 0} 
and the eigenvalue with multiplicity n, since —Q22 has its spectrum in {Re A > 0}. The 
classical perturbation theory as in [MaZlj shows that for rj~2 small, rp A has n eigenvalues 
of order rj~2, close to the spectrum of A^ with error 0{r]~ l ). Thus, for 7/ > large, 772 A 
has r + dim S{A\ 1 ) eigenvalue in {Re A < 0}, proving (I9.17|) . 

The proof of the Proposition is now complete. □ 

9.3 Finite dimensional approximations 

To treat the infinite-dimensional case, we proceed by finite-dimensional approximations. 
Let Q = Qm_ and K_ = Km_ denote the operators Qu and K\j evaluated at the equilibrium 
M^M(u), so that 



(9.19) A=(A)*, Q = (Q)* 




Q 



.22 



with 

(9.20) Q 99 > eld, c > 0. 
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Ku K12 
K21 



Moreover, the Kawashima multiplier has the form 

(9.21) K = -(K)* = d 

Thanks to (|9.20p the condition (|4.14p is satisfied for 9 small enough as soon as 

(9.22) Re (K n A n + K 12 A 21 ) > eld, c> 0. 
Consider an increasing sequence of finite dimensional subspaces 

(9.23) V r C V r +i, UV r dense in V. 

Similarly, let HL, = Uffi V r . Let U r (U) denote the orthogonal projector onto EL., so that 

(9.24) n r = n*, 
and define 

(9.25) A r (U) = U r (U)AU r (U), Q r (U) = U r (U)QuU r (U). 

Lemma 9.5. (i) A r (U) is symmetric. 

(ii) Q r {U) is uniformly definite negative on {0} © V r . 
(Hi) With K r : U r (U)KjjIl r (U) , the Kawashima condition 

(9.26) Re K r A r - Q r > 7H 

is uniformly satisfied for U in a neighborhood of M_ and r large enough. 

Proof, (i) follows by (|9.24|) and symmetry of A on H. 

On {0} © V r , Q r (U) is a perturbation of n r Q ir r where 7r r is the (usual) orthogonal 
projection onto V r , with 

R r ( u ) = ( u 

To prove (I9.26p . is is sufficient to prove the property at U = M_. Restricting to V r by 
7r r , one has 

An Au^r \ TS — a ( —11 K-U^r 



^ V \KrA.21 7T rA22 7T r) ' ^' ^ \ 7r rK.2l 

Note that 

I£llAll + K_l2 7T rA2i 

is an n-dimensional perturbation of K_nAii + ^12^.21 whose real part is definite positive. 
Therefore, for r large enough, 

Re (K r A r ) u = Re (K U A U + K 12 TT r A 21 ) > eld 

with c independent of r. Since 

TTrQ 22 TTr > Clld On V r 

uniformly in r, and since the other blocks of K r A r are uniformly 0(9), the condition (|9.26l) 
is satisfied for r large enough and 9 sufficiently small. □ 
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Corollary 9.6. On H r the equation 

(9.27) A r (u NS )d x U - Q r {u NS )U = (Q , I ■ u(0) = 

is well posed, and there are uniform estimates in r, for r sufficiently large. 
Corollary 9.7. On M the equation 

(9.28) Ad x U - Q(u NS )U = r\ I ■ u(0) = 
is well posed. 

9.4 Proof of Proposition 16.21 

Let (Ct' 11 )^ denote the inverse operator of C% ,n denned by (|9.ip . for sufficiently small rj > 0. 
The uniform bound (|9.3h . and weak compactness of the unit ball in H 2 , for F 6 y, we 
obtain existence of a weak solution U G H 2 of 

(9.29) ClU = F := Q , £ e ■ u(0) = 0, 

along some weakly convergent subsequence. Proposition 18.21 implies uniqueness in H 2 for 
this problem, therefore the full family converges, giving sense to the definition 

(9.30) (£ £ jt = l im (£^)t 

77 — i-0 

acting from y to H 2 . 

For F G S £ s, the uniform bounds ()9.3p imply that the limit (C%yU S il| s and satisfies 

same estimate. By density, the operator {C%)^ extends to / G -fff^ 1 and gf G with 
(£ £ Jt FGjff s 5 . 

The sharp bound (|6.14p and (|6.15p now follow immediately from Propositions 18.21 and 
The proof of Proposition 16.21 is now complete. 



10 Other norms 

We now briefly discuss the modifications needed to obtain the full result of Theorem 11.11 
10.1 Pointwise velocity estimates 

Proof of Theorem \l.l\ (H 1 / 2 ). To obtain pointwise bounds with respect to velocity, we carry 
out the same argument as in the proof of Proposition 15.41 substituting in place of the L 2 
norm | • | in £, the weighted H s (Sobolev) norm 

s 

\f\ s :=J2C' k \dlf\ 2 , 

k=0 
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C > sufficiently large, similarly as we did for the x-variable in order to get pointwise 
bounds in x. 

We have only to observe that differentiating the linearized equations in £ gives the same 
principal part applied to the ^-derivative of U, plus commuator terms. Since commutator 
terms, both for the linearized collision operator L and the transport operator A are of one 
lower derivative in £ and also one lower factor in (£) (straightforward computation differ- 
entiating |£ — £1, respectively) for the hard-sphere case, we easily find that commutator 
terms are absorbable for C > sufficiently large by lower order estimates already carried 
out. 

Thus, we obtain all the same estimates as before and the argument closes to give the 
same result in the stronger norm | • | s . (Note: a detail is to observe that truncation errors of 
the approximate solution are of the same order in the Sobolev norm, which follows by the 
corresponding property of the Maxwellian.) Applying the Sobolev embedding estimate in 
£, we obtain (j!.17[) for r) = 1/2, which evidently implies the same estimate for 77 > 1/2. □ 

10.2 Higher weights 

Proof of Theorem li.il (H 3 ). To extend our results from EI 2 to HP, we use a simple bootstrap 
argument together with the key observation that the HP norm of Pu/ is controlled (by 
equivalence of finite-dimensional norms) by the H2 estimates already obtained. Namely, 
starting similarly as in (|4.24p with the equation Ad x — L a U = F, F^F = f, FyF = g, 
we find, taking the HP-inner product of U against this equation and applying the result of 
Proposition 13.61 and recalling that A is formally self-adjoint in HP, we obtain the estimate 



Differentiating the equations k times and taking the inner product with d^U, we find, 
similarly, the higher-derivative estimate 



Specializing now to the case (|6.8p . (|6.9p . and bounding the HP norm of Pull by a 
constant times the Ha bound obtained already in our previous analysis, we recover the key 
bounds (|6.14p - (|6.15p of Proposition HOI in the general space M s . With this bound, the entire 
contraction mapping argument goes through in HP, since this relies only on boundedness 
estimates on A, Q already obtained, the estimate (|5.2p (still valid in HP), and the linearized 
estimates Q6.14j> — f]6.15 j> . yielding (|1.17p (i) and (ii) as claimed, for any 77 > 0. 

The estimate ()1.17p (iii) then follows by Remark 13 . 71 estimating decay in velocity £ of the 
approximating profile JnS- D 

Remark 10.1. We emphasize that L is not approximately self-adjoint with respect to 
HP, s » 1/2, and, likewise, the splitting HP = U © V s using projectors P[/ and Fy is not 
orthogonal in this norm. For this reason, we obtain term ||PuC/||^fc in the righthand side 



(10.1) 




(10.2) 



\\VyU\\ H k < C{\\f\\ H k + \\g\\ H , + ||Puf/|UO- 

e,o £,5 £,5 e,6' 
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of (|10.2p and not e||PuC/|| H k as in the H 1 / 2 computations. The missing e factor was crucial 

in closing the argument in H 1 / 2 and estimating Put/. However, with Put/ already bounded 
it is no longer needed, since our final estimates make no distinction between Put/ and Pyt/ 
components; that is, the lost e factor is needed only to close the loop between microscopic 
and macroscopic estimates, and not to bound FyU in terms of Put/. 

11 Other potentials 

Finally, we briefly indicate the changes needed to accomodate general hard cutoff potentials. 
Recall [CNl IGPS} that these give structure L = -v(£) + /C, where v ~ (£)P, < /3 < 1, 
and K, is compact, and similarly for Q. Dividing by v ~ as before, we can thus obtain 
Q, L bounded, but this leaves A unbounded. Nonetheless, a closer look shows that the 
Kawashima compensator K as constructed is still bounded, the key point. For, examining 
A12, we see that it decays as a polynomial in £ times full Maxwellian, so is clearly bounded 
in HP for s < 1. 

Since the norm of A does not enter except through the good term KA, our basic micro- 
estimates therefore still survive. Of course, the macro-estimates, since finite-dimensional, 
survive as well. (This follows by the same estimate that shows that K as constructed is 
bounded; that is, one has only to check that A\2 and A21 entries remain bounded, thanks 
to Maxwellian rate decay.) Thus, the argument goes through as before, also for these more 
general potentials. 
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